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ABSTRACT 

The  nearby  star  a  Oph  (Ras  Alhague)  is  a  rapidly  rotating  A5IV  star  spinning  at  ~  89%  of  its  breakup  velocity. 
This  system  has  been  imaged  extensively  by  interferometric  techniques,  giving  a  precise  geometric  model  of  the 
star’s  oblateness  and  the  resulting  temperature  variation  on  the  stellar  surface.  Fortuitously,  a  Oph  has  a  previously 
known  stellar  companion,  and  characterization  of  the  orbit  provides  an  independent,  dynamically  based  check  of 
both  the  host  star  and  the  companion  mass.  Such  measurements  are  crucial  to  constrain  models  of  such  rapidly 
rotating  stars.  In  this  study,  we  combine  eight  years  of  adaptive  optics  imaging  data  from  the  Palomar,  AEOS, 
and  CFHT  telescopes  to  derive  an  improved,  astrometric  characterization  of  the  companion  orbit.  We  also  use 
photometry  from  these  observations  to  derive  a  model-based  estimate  of  the  companion  mass.  A  fit  was  performed 
on  the  photocenter  motion  of  this  system  to  extract  a  component  mass  ratio.  We  find  masses  of  2.40+ q2337  M©  and 
0.85+qq6  Mq  for  a  Oph  A  and  a  Oph  B,  respectively.  Previous  orbital  studies  of  this  system  found  a  mass  too 
high  for  this  system,  inconsistent  with  stellar  evolutionary  calculations.  Our  measurements  of  the  host  star  mass 
are  more  consistent  with  these  evolutionary  calculations,  but  with  slightly  higher  uncertainties.  In  addition  to  the 
dynamically  derived  masses,  we  use  IJHK  photometry  to  derive  a  model-based  mass  for  a  Oph  B,  of  0.77  zb  0.05 
Mq  marginally  consistent  with  the  dynamical  masses  derived  from  our  orbit.  Our  model  fits  predict  a  periastron 
passage  on  2012  April  19,  with  the  two  components  having  a  50  mas  separation  from  2012  March  to  May.  A  modest 
amount  of  interferometric  and  radial  velocity  data  during  this  period  could  provide  a  mass  determination  of  this 
star  at  the  few  percent  level. 

Key  words:  instrumentation:  adaptive  optics  -  methods:  data  analysis  -  stars:  individual  (HIP86032)  -  techniques: 
image  processing 


1.  INTRODUCTION 

Although  the  binary  properties  of  solar-type  stars 
(Duquennoy  &  Mayor  1991;  Metchev  &  Hillenbrand  2009)  and 
lower  mass  late-type  stars  (Fischer  &  Marcy  1992;  Reid  &  Gizis 
1997)  are  becoming  increasingly  well  understood,  studies  of  the 
higher  end  of  the  mass  spectrum  are  crucial  for  a  broader  under¬ 
standing  of  stellar  multiplicity.  Indeed,  recent  simulations  sug¬ 
gest  that  A  stars  may  possess  a  greater  abundance  of  both  stellar 
and  planetary  mass  companions  (Kratter  et  al.  2010,  J.  Crepp 
2010,  private  communication).  However,  the  multiplicity  char¬ 
acteristics  of  these  more  massive  stars  have  only  recently  started 
to  be  surveyed  (Shatsky  &  Tokovinin  2002;  Kouwenhoven  et  al. 
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2005),  and  more  extensive  surveys  are  underway  to  place  limits 
on  the  nature  of  massive-star  multiplicity.  In  recent  years,  high 
contrast  imaging  techniques  (Oppenheimer  &  Hinkley  2009; 
Absil  &  Mawet  2010)  such  as  adaptive  optics  (hereafter  “AO”) 
have  matured  significantly,  allowing  these  studies  to  proceed. 
The  nearby  (14.68  pc)  A5IV  star  a  Oph  (Ras  Alhague)  has 
a  well-known  companion  (Wagman  1946;  Lippincott  &  Wag- 
man  1966;  Gatewood  2005)  with  a  8.62  yr  period,  well  estab¬ 
lished  over  several  decades  of  monitoring  and  first  resolved  by 
McCarthy  (1983).  But  a  fuller  characterization  of  the  compan¬ 
ion  has  not  been  possible,  and,  aside  from  the  period,  the  orbital 
parameters  have  been  loosely  constrained  (Kamper  et  al.  1989; 
Augensen  &  Heintz  1992;  Gatewood  2005). 

Characterization  of  the  binary  nature  of  A-star  systems  be¬ 
comes  doubly  interesting  when  the  host  stars  are  also  rapidly 
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Table  1 

Table  of  Astrometry 


Hinkle y  et  al. 


Epoch 

MJD 

r  (mas) 

Or 

RA. 

°P.A. 

x,  y  (mas) 

Band 

Telescope  (Instrument) 

Ref 

1 

2445157.9 

— 380y  ±  50 

K,  L 

3.8m  KPNO  May  all 

1 

2 

2445244.8 

—65 Ox  ±  50 

K,  L 

3.8m  KPNO  Mayall 

1 

3 

2445424.2 

— 340y  ±  50 

K,  L 

3.8m  KPNO  Mayal 

1 

4 

2445804.0 

— 430y  ±  50 

K,  L 

3.8m  KPNO  Mayall 

1 

5 

2447045.8 

— 400x  ±  50 

K,  L 

3.8m  KPNO  Mayall 

1 

6 

2451304.5 

770 

±40 

243.7 

±3?0 

2.15  /xm 

3.6m  La  Silla  (SHARP  II) 

2 

7 

2452473.8 

470 

±20 

233.2 

±2?4 

/ 

AEOS  (Vislm) 

3 

8 

2452760.0 

<275 

±3 

n/a 

/ 

AEOS  (Vislm) 

3,4 

9 

2453168.5 

303 

±33 

253.0 

±6?3 

H 

AEOS  (Lyot  Project) 

5 

10 

2454253.9 

776.5 

±2.1 

244.6 

±0?4 

J,  H,  K 

Palomar  (PHARO) 

6 

11 

2454263.5 

765 

±20 

243.7 

±  1?5 

H 

AEOS  (Lyot  Project) 

5 

12 

2454637.5 

787.8 

±2.8 

240.6 

±0?4 

K,Bry 

Palomar  (PHARO) 

6 

13 

2454657.5 

790 

±20 

239.5 

±  1?4 

J,H 

Palomar  (Project  1640) 

7 

14 

2454963.8 

756.6 

±7.5 

239.3 

±  1?2 

Pa^  (1.28  /xm) 

CFHT  (PUEO) 

8 

15 

2455002.5 

748.5 

±2.7 

238.3 

±0?4 

K 

Palomar  (PHARO) 

6 

Notes.  (1)  McCarthy  (1983)  and  Kamper  et  al.  (1989)  provide  one-dimensional  speckle  measurements,  (2)  Boccaletti  et  al.  (2001),  (3)  Roberts  &  Neyman 
(2002),  (4)  Upper  limit  to  orbit,  (5)  Sivaramakrishnan  et  al.  (2007),  (6)  Hayward  et  al.  (2001),  (7)  Hinkley  et  al.  (2008,  201 1),  (8)  Rigaut  et  al.  (1998). 


rotating.  Recent  interferometric  imaging  of  several  rapidly 
rotating  A  stars  have  revealed  imaging  of  their  surfaces 
(Monnier  et  al.  2007;  Zhao  et  al.  2009).  a  Oph,  a  prototypical 
rapid  rotator,  is  rotating  at  ~89%  of  its  breakup  velocity.  In¬ 
terferometric  surface  imaging  of  this  star  clearly  shows  its  pro¬ 
nounced  oblateness:  the  star’s  equatorial  radius  is  20%  larger 
than  its  polar  radius,  corresponding  to  a  ~1840  K  tempera¬ 
ture  differential  (Zhao  et  al.  2009).  This  2-3  M0  star  has  also 
been  the  target  of  extensive  asteroseismic  monitoring  carried 
out  by  the  MOST  satellite,  discovering  ~50  pulsational  modes 
(Monnier  et  al.  2010). 

Derivation  of  the  component  masses  for  this  system  is  a  crit¬ 
ical  check  against  investigations  of  the  host  star’s  rapid  rotation 
and  asteroseismology.  Discrepancies  between  the  dynamically 
derived  mass  of  the  host  star  and  that  allowed  by  models  of 
rapid  rotators  may  illuminate  potential  complications  with  rota¬ 
tor  models,  e.g.,  those  that  do  not  take  into  account  differential 
rotation  of  the  star.  For  all  these  reasons,  deriving  a  more  ac¬ 
curate  orbit  is  crucial  for  future  studies  to  reconcile  the  orbital 
dynamics  of  the  companion  star  with  the  physics  dictating  the 
rotation  of  the  host  star.  In  this  paper,  we  assemble  data  from  a 
number  of  high-contrast  imaging  programs  over  the  past  eight 
years.  We  take  advantage  of  these  recent  measurements  along 
with  archival  data  to  further  constrain  the  orbit  of  the  companion 
as  well  as  the  system  mass. 

2.  ASTROMETRIC  OBSERVATIONS 

Table  1  presents  astrometric  measurements  of  a  Oph 
obtained  using  several  instruments.  Epochs  1  through  5 
are  previously  measured  one-dimensional  speckle  data  from 
Kamper  et  al.  (1989).  For  the  sixth  epoch  of  observations,  we 
adopt  the  radial  separation  (770±  40  mas)  and  position  angle 
(243?7  zb  3?0)  measurements  reported  in  Boccaletti  et  al.  (2001). 
Boccaletti  et  al.  (2001)  mistakenly  reported  the  position  angle 
in  degrees  measured  west  of  north.  Gatewood  (2005)  noticed 
the  problem,  but  mistakenly  assumed  it  was  a  quadrant  prob¬ 
lem  and  changed  the  value  by  180°.  We  use  the  correct  mea¬ 
surement  position  angle  of  243°  (A.  Boccaletti  2010,  private 
communication). 

Epochs  7  and  8  are  comprised  of  visible  AO  measure¬ 
ments  obtained  from  the  3.63  m  AEOS  telescope  on  Haleakala, 


Figure  1.  Visual  orbit  of  the  a  Oph  system.  The  solid  line  is  the  best-fit  solution 
while  the  gray  band  shows  the  allowed  orbits  based  on  Monte  Carlo  sampling  of 
the  parameter  uncertainties.  The  primary  is  located  at  (0,0)  and  the  astrometric 
data  from  Table  1  are  plotted  with  the  positional  uncertainties.  East  is  left  and 
north  is  up. 


Hawaii  (Roberts  &  Neyman  2002;  Roberts  et  al.  2005). 
Epochs  9  and  11  were  obtained  using  “The  Lyot  Project” 
(Sivaramakrishnan  et  al.  2007;  Hinkley  et  al.  2007,  2009; 
Leconte  et  al.  2010),  a  diffraction-limited  classical  Lyot  coro- 
nagraph  (Lyot  1939;  Sivaramakrishnan  et  al.  2001)  working  in 
the  infrared  and  recently  decommissioned  at  AEOS. 

The  thirteenth  epoch  of  observations  of  the  a  Oph  system  was 
obtained  using  a  recently  commissioned  coronagraph  integrated 
with  an  integral  field  spectrograph  (IFS)  spanning  the  J  and  H 
bands  (1 .06  fim-l  .76  /xm)  on  the  200  in  Hale  Telescope  at  Palo- 
mar  Observatory  (Hinkley  et  al.  2008,  2011).  This  instrument 
package,  called  “Project  1640”  is  mounted  on  the  Palomar  AO 
system  (Dekany  et  al.  1998;  Troy  et  al.  2000)  and  is  a  dedicated 
high  contrast  imaging  instrument  providing  low  resolution  spec¬ 
tra  30-60  with  a  lenslet-based  IFS  and  Apodized-Pupil 

Lyot  coronagraph  (Soummer  2005),  an  improvement  of  the  clas¬ 
sical  Lyot  coronagraph  (Sivaramakrishnan  et  al.  2001).  Epoch 
14  observations  were  made  at  the  Canada-France-Hawaii 
Telescope  (Rigaut  et  al.  1998).  The  observations  comprising 
epochs  10,  12,  and  15  were  obtained  with  the  Palomar  AO  sys¬ 
tem  and  the  PHARO  infrared  camera  (Hayward  et  al.  2001). 
Astrometric  points  from  epochs  6  through  15  are  shown  in 
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arcsec  arcsec  arcsec 

Figure  2.  Examples  of  data  taken  at  the  three  epochs  listed  at  the  upper  left  hand  corner  of  each  panel.  The  left  and  middle  images  were  taken  with  a  near-IR  camera 
with  a  coronagraph  at  the  AEOS  telescope.  The  right  hand  panel  is  a  single  slice  extracted  from  a  data  cube  from  the  high  contrast  integral  field  spectrograph  at 
Palomar  Observatory  (Hinkley  et  al.  2008,  201 1).  The  companion,  a  Oph  B,  is  marked  with  a  white  arrow  in  all  three  epochs,  and  is  shown  near  apastron  in  the  third 
panel. 


Table  2 

Binary  Orbit  and  Derived  Properties 

Parameter 

Value 

From  visual  orbit  only 

Semi-major  axis  a  (mas) 

Eccentricity 

Inclination  (deg) 
co  (deg) 

Q  (deg) 

Period  (days) 

Epoch  of  periastron  (JD) 

427^ 

0.92  ±  0.03 
125+| 

162  ±14 

232  ±9 
3148.4  (fixed) 
2452888  ±  53 

From  astrometric  photocenter  motion 
(using  visual  orbit  above) 

Mass  ratio  (primary/secondary) 

Parallax  (mas) 

Fixed  proper  motion:  E,  N  (mas  yr-1) 

System  mass  (M0) 
a  Oph  A  (Af0) 
a  Oph  B  (M0) 

2.76  +423 
69.1  (fixed) 
123.3,  -227.0 
3  25+-35 

2.40+f7 

0.85+“ 

Figure  1,  and  representative  examples  of  data  from  epochs  9, 
11,  and  13  are  shown  in  Figure  2. 

3.  ORBITAL  CHARACTERIZATION 

We  fitted  an  orbit  to  the  astrometric  data  in  Table  1  includ¬ 
ing  the  speckle  data  from  Kamper  et  al.  (1989). 19  We  used 
both  a  Levenberg-Marquardt  method  (using  IDL  procedure 
MPFIT  developed  by  C.  Markwardt20),  and  also  the  IDL  rou¬ 
tine  AMOEBA  (Press  et  al.  1992),  arriving  at  the  same  global 
solution.  Comparing  our  predicted  astrometry  with  the  data,  we 
find  a  reduced  x2  °f  0-7.  The  best- fit  orbital  elements  can  be 
found  in  Table  2,  where  we  have  fixed  the  quadrant  of  co  using 
the  radial  velocity  data  of  Kamper  et  al.  (1989).  In  order  to  cal¬ 
culate  errors  for  the  orbital  elements,  we  carried  out  a  Monte 
Carlo  simulation,  allowing  all  parameters  to  vary  in  the  fitting 

19  We  needed  to  flip  the  signs  on  x,  y  from  this  work  to  be  consistent  with 
other  astrometry  measurements. 

20  http://cow.physics.wisc.edu/~craigm/idl/idl.html 
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Figure  3.  Relative  photocenter  positions  of  a  Oph  from  MAP  observations 
(Gatewood  1987)  reported  by  Gatewood  (2005)  after  subtracting  out  the  system 
proper  motion  and  an  estimate  of  the  center  of  mass.  With  the  orbit  fit  shown  in 
Figure  1,  the  component  mass  ratio  is  obtained  by  fitting  the  photocenter  motion 
shown  here. 


except  for  the  period,  which  we  fixed  to  8.62  yr.  The  results  of 
the  error  analysis  can  be  found  in  Table  2. 

Using  the  full  co  variant  set  of  parameters  from  the  Monte 
Carlo  study,  we  were  then  able  to  determine  the  component 
mass  ratio  by  fitting  the  photocenter  motion  measured  with  the 
MAP  instrument  (Gatewood  2005).  In  order  to  carry  out  a  fit  to 
the  absolute  position  of  the  photocenter,  we  needed  to  constrain 
a  few  important  parameters  using  independent  sources.  First, 
we  adopted  a  proper  motion  of  123.3  mas  yr-1  (east)  and 
—227.0  mas  yr-1  (north)  from  the  FK5  catalog.  Note  that 
other  proper  motion  catalogs  with  shorter  time  baselines  contain 
significantly  different  values  for  a  Oph,  presumably  having  been 
corrupted  by  orbital  motion  of  the  primary.  We  also  adopted  a 
parallax  of  69.1  mas  from  Martin  &  Mignard  (1998).  Note 
that  the  errors  in  these  three  quantities  have  been  neglected  in 
the  following  analysis  (errors  from  our  visual  orbit  dominate 
uncertainties  in  the  mass  ratio).  Figure  3  shows  these  relative 
photocenter  offsets  after  subtracting  out  the  system  proper 
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Table  3 

Apparent  Magnitudes  for  the  a  Oph  System 


Member 

/ 

J 

H 

K 

a  Oph  Aa 

1.90 

1.83  ±  0.28 

1.72  ±  0.18 

1.68  ±  0.21 

a  Oph  Bb 

6.72  ±  0.10 

6.02  ±  0.31 

5.44  ±  0.21 

5.25  ±  0.24 

Notes. 

a  I- band  photometry:  Monet  et  al.  (2003),  JHK :  Cutri  et  al.  (2003) 
b  /-band  photometry  was  obtained  using  the  AEOS  Vislm  instrument  and  JHK 
photometry  with  the  PHARO  infrared  camera  at  Palomar. 

motion  and  an  estimate  of  the  center  of  mass.  Lastly,  we  had 
to  assume  a  flux  ratio  at  the  effective  wavelength  of  the  MAP 
experiment  Gatewood  (1987),  similar  to  an  R-band  filter.  Based 
on  multi-wavelength  detections  of  the  companion  at  /,  /,  H ,  and 
K  bands,  we  have  estimated  an  R- band  flux  ratio  of  5.3  mag 
(factor  of  130)  assuming  the  companion  is  a  mid/late  K  star 
(See  Section  4.1).  Using  the  results  from  the  two  analyses,  we 
constrain  the  masses  to  be  2.40+q  37  Mq  and  0.85/°0°064  M0  for 
a  Oph  A  and  a  Oph  B,  respectively.  Table  2  contains  the  results 
of  this  orbital  fitting,  along  with  the  full  set  of  the  fit  parameters. 

4.  ANALYSIS 

4.1.  Companion  Spectral  Energy  Distribution 

Table  3  presents  the  broadband  apparent  magnitudes  for  a 
Oph  B  obtained  with  the  PHARO  infrared  camera  at  Palo¬ 
mar,  as  well  as  the  AEOS  Vislm  camera.  The  spectral  energy 
distribution  for  a  Oph  B  is  shown  in  Figure  4,  accounting 
for  the  system’s  0.834  zb  0.024  distance  modulus  (Gatewood 
2005).  The  broadband  magnitudes  listed  in  Table  3  are  the 
result  of  a  simultaneous  photometric  fit  to  each  binary  com¬ 
ponent  which  returns  position  angle,  relative  separation,  as 
well  as  relative  brightnesses.  Figure  4  shows  /-,  /-,  //-,  BrK-, 
and  K- band  photometry.  Also  shown  are  broadband  J-  and 
H- band  photometric  values  obtained  with  the  Project  1640  IFS 
at  Palomar  Observatory.  The  Project  1640  observations  were 
calibrated  using  reference  star  spectra  obtained  from  the  IRTF 
spectral  library  (Cushing  et  al.  2005;  Rayner  et  al.  2009)  in  the 
same  manner  as  described  in  Hinkley  et  al.  (2010).  However, 
these  observations  were  gathered  during  the  first  light  observa¬ 
tions  with  this  new  instrument,  and  hence  the  calibration  data 
set  was  obtained  before  a  full  understanding  of  an  effective 
calibration  process  was  mature. 

Also  shown  in  Figure  4  are  theoretical  broadband 
and  K- band  fluxes  based  on  the  models  of  Siess  et  al.  (2000).  The 
curves  in  Figure  4  are  interpolations  between  these  values.  At 
a  spectral  resolution  of  30-60,  late-G  to  early-M  star  spectra 
are  very  similar,  preventing  robust  discrimination  between 
them.  Rather  than  plotting  template  spectra  over  the  data  for 
comparison,  we  have  instead  chosen  to  use  these  model-based 
interpolated  mass-luminosity  curves  to  allow  discrimination  of 
the  best- fit  mass  value  for  a  Oph  B.  The  best- fit  curve  to  the 
measured  fluxes  corresponds  to  a  mass  of  0.77  Mq.  Taking  into 
account  our  photometric  errors  as  well  as  the  0.024  mag  distance 
modulus  uncertainty  quoted  by  Gatewood  (2005),  we  deduce  an 
uncertainty  of  ±0.05  M0  in  our  best-fit  value.  This  value  of  0.77 
±  0.05  Mq  is  very  similar  to  the  0.78  ±  0.058  M0  reported 
by  Gatewood  (2005),  and  within  1.3a  of  our  0.85+J 04  Mq 
dynamical  mass  discussed  in  Section  3.  Moreover,  in  this  mass 
regime  (~0.7-0.8  Mq),  Hillenbrand  &  White  (2004)  have 
demonstrated  that  mass-luminosity  models,  including  those 
produced  by  Siess  et  al.  (2000),  are  only  inaccurate  at  the  few 


Figure  4.  Broadband  photometry  for  a  Oph  B  from  three  of  the  observing 
programs  described  in  this  paper.  Also  shown  are  three  theoretical  broadband 
/-,  J-,  H-,  and  K- band  flux  curves  based  on  the  mass-luminosity  models  of  Siess 
et  al.  (2000).  Each  curve  is  an  interpolation  between  these  values.  The  best 
fitting  model  is  shown  (solid  line — 0.77  Mq),  as  is  the  dynamically  determined 
mass  of  0.85  Mq,  and  a  0.65  Mq  curve  for  reference.  The  horizontal  error  bars 
reflect  the  width  of  each  bandpass. 


percent  level.  Due  to  the  aforementioned  uncertainties  in  the 
calibration  of  the  Palomar  Project  1640  photometry,  only  the 
Palomar  PHARO  data  and  the  AEOS  I- band  data  were  used 
in  the  fit.  Figure  4  also  shows  a  curve  corresponding  to  the 
dynamical  best- fit  mass  value  of  0.85  Mq,  as  well  as  a  curve 
corresponding  to  0.65  Mq  for  reference.  Comparing  our  derived 
near-infrared  ( J  —  H )  and  ( H  —  K )  colors  with  those  documented 
in  Bessell  &  Brett  (1988)  put  the  spectral  type  of  the  companion 
between  a  K5V  and  a  K7V. 

4.2.  Spin-Orbit  Mutual  Inclinations 

Our  measurement  of  the  orbital  plane  of  the  a  Oph  binary 
system  can  be  compared  to  the  angular  momentum  of  the  rapidly 
rotating  primary  star.  This  is  analogous  in  practice  to  measuring 
the  relative  orbital  inclinations  of  triple  or  quadruple  system 
(e.g.,  Muterspaugh  et  al.  2008)  which  can  reveal  hints  to  the 
formation  and/or  evolution  of  system  dynamics  (e.g.,  Sterzik 
&  Tokovinin  2002).  In  order  to  carry  out  this  comparison,  we 
convert  the  spin  parameters  from  Zhao  et  al.  (2009)  into  an 
equivalent  Cl  =  216?  12  ±  1?23,  i  =  92?30  ±  0?43  that  would 
characterize  a  circular  orbit  in  the  star’s  equatorial  plane.  Note 
that  these  values  are  affected  by  orbital  degeneracies  since  we 
have  not  incorporated  any  surface  or  orbital  velocities  into  the 
estimate.  However,  because  the  primary  spin  geometry  is  nearly 
edge-on,  these  degeneracies  do  not  matter  much.  We  find  that  the 
mutual  angle  of  inclination  is  O  =  36°  ±  6°  (or  O  =  144°  ±  6° 
for  the  opposite  spin  polarity).  Interestingly,  this  lies  close  to  the 
predicted  peak  inclination  (40°)  that  arises  from  Kozai  cycles 
in  the  context  of  an  inner  short-period  binary  (Fabrycky  & 
Tremaine  2007),  echoing  findings  of  Muterspaugh  et  al.  (2008). 
However,  the  primary  a  Oph  A  has  no  known  inner  companion 
and  so  this  mechanism  does  not  seem  to  be  applicable. 

5.  DISCUSSION  AND  CONCLUSIONS 

We  have  obtained  an  orbital  fit  for  the  a  Oph  system  and 
derived  component  masses  of  2.40+-2337  and  0.85t'Q64  M©  for 
a  Oph  A  and  a  Oph  B,  respectively.  Our  estimation  of  the 
mass  of  a  Oph  A  is  lower  than  the  2.84  ±  0.19  Mq  value 
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quoted  by  Gatewood  (2005).  However,  our  Monte  Carlo  ap¬ 
proach  yields  more  conservative  estimates  of  the  uncertainties, 
overlapping  with  other  works.  Specifically,  the  2. 1  M©  value  for 
a  Oph  A  quoted  by  Zhao  et  al.  (2009)  and  measured  using  an 
H-R  diagram  overlaid  with  isochrones  is  allowed  at  this  work’s 
lcr  level,  while  it  is  only  allowed  at  the  3.5cr  level  in  the  Gate- 
wood  (2005)  study.  Our  work  demonstrates  that,  especially  with 
a  high  eccentricity  orbit,  careful  estimates  of  stellar  parameters 
are  essential. 

We  anticipate  that  this  dynamical  estimation  of  the  mass  of 
a  Oph  A  will  aid  forthcoming  studies  focused  not  only  on 
modeling  the  rapid  rotation  of  a  Oph  A,  but  also  of  the  host 
star’s  asteroseismological  properties  (Monnier  et  al.  2010).  The 
mass  value  derived  here  may  help  to  distinguish  between  models 
that  assume  solid-body  rotation  of  the  star  or  have  some  degree 
of  differential  rotation  beneath  the  photosphere. 

Moreover,  we  have  compiled  the  most  complete  set  of 
photometric  data  yet  for  a  Oph  B  from  several  observing 
programs  to  obtain  an  estimate  of  the  companion  mass  through 
the  use  of  theoretical  models.  The  0.77  zb  0.05  M©  estimation 
of  the  companion  mass  derived  photometrically  is  marginally 
consistent  with  the  0.85+- q64  M©  value  for  a  Oph  B.  Our  infrared 
colors  and  the  corresponding  best-fit  mass-luminosity  curve 
are  consistent  with  a  mid/late-K  dwarf  classification  for  the 
companion. 

Finally,  based  on  our  fit,  a  Oph  B  will  be  passing  through 
periastron  on  2012  April  19,  with  an  uncertainty  on  this  date 
of  ±  53  days.  However,  the  separation  between  the  two  objects 
should  be  ~50  mas  from  2012  March  24  to  2012  May  21. 
Such  a  separation  will  be  ideal  for  existing  state-of-the-art 
interferometers  such  as  the  CHARA  Array  or  the  Very  Large 
Telescope  Interferometer.  Intensive  astrometric  monitoring  of 
this  system  and  new  radial  velocity  observations,  especially 
during  periastron,  can  help  to  constrain  the  mass  of  a  Oph  A  to 
within  a  few  percent. 
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